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Abstract 

Dimensional reduction and effective field theory methods are applied to the Two Higgs Doublet 
Model at finite temperature. A sequence of two effective three-dimensional field theories which are 
valid on successively longer distance scales is constructed. The resulting Lagrangian can be used to 
study different aspects of the phase transition in this model as well as the sphaleron rate immediately 
after the phase transition. 
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1 Introduction 



The electroweak phase transition (EWPT) has been the subject of intense investigation in recent 
years, largely due to its possible role in generating the baryon asymmetry of the present Universe [1- 
2] (see also for a detailed review). If the electroweak phase transition is of first order, it 
proceeds through bubble nucleation, and the baryon asymmetry is produced as the bubbles expand 
and the Universe is far from equilibrium [1-2]. Moreover, if the baryon asymmetry produced 
during the phase transition has survived until today, baryon number violating processes (sphaleron 
processes) must have been suppressed immediately after the phase transition §1. This requires 
that the electroweak phase transition is strongly first order [|J. It is now well established that this 
requirement is not met in the Standard Model (SM); for realistic values of the Higgs mass, the 
phase transition is either too weakly first order to suppress the sphaleron processes, or it is second 
order, or there is no phase transition at all |J. Furthermore, it is not clear that the amount of 
CP violation in the SM is sufficient. 

The fact that baryogensis is ruled out in the Standard Model suggests the investigation of 
extensions of the Standard Model such as the Minimal Supersymmetric Standard Model (MSSM) 
and the Two Higgs Doublet Model (2HDM) [[]]. Both these theories have additional sources of 
CP violation. The main objective is to find regions in parameter space where the phase transition 
is strongly first order, so that the excess of baryons produced during the phase transition is not 
washed out by sphaleron processes immediately after the phase transition. 

The electroweak phase transition in the Standard Model [7-9] as well as in the Two Higgs 
Doublet Model [10-14] has been studied using resummed perturbation theory. The strength of the 
EWPT in the Standard Model weakens as the Higgs mass increases [7-9]. However, the resummed 
loop expansion breaks down for large (realistic) Higgs masses (see e.g. Ref. || for a discussion of 
the validity of the resummed perturbation expansion), and so one must employ nonperturbative 
methods in order to discriminate between a weakly first order and a second order phase transition. 

The electroweak phase transition in the SM has also been investigated by lattice simulations di- 
rectly in four dimensions [15-17], renormalization group techniques [18-19] and the e-expansion p0| . 



These methods yield results for the quantities characterizing the phase transition that are in qual- 
itative agreement with the perturbative treatment. 

Significant progress in the study of phase transitions has been made by applying the methods 



of dimensional reduction [21-24] and effective field theory [25]. The idea is to integrate out the 
nonzero bosonic modes as well as the fermionic modes which decouple from the static modes at high 
temperature [21-24]. One is then left with an effective three-dimensional field theory of the zero 
modes. Since the nonstatic modes have masses of order T, the process of dimensional reduction 
is free of infrared problems, and in weakly coupled theories this can normally be carried out 
perturbatively. In nonabelian gauge theories the effective three-dimensional theory contains two 



momentum scales [BBJ . The scale gT is provided by the temporal component of the gauge field, and 
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the scale g 2 T is provided by the spatial components of the gauge field. Moreover, in theories with 
a single Higgs multiplet, the scalar mass is normally of order g 2 T for temperatures close to T c . For 
theories with more than one Higgs multiplet, the masses of the additional scalar fields are generally 
of order gT. In either case it proves useful to construct a second effective field theory by integrating 
out the timelike component of the gauge field (and possibly some scalar fields) as shown by Farakos 



et al. |27| and by Braaten and Nieto p8fl . This approach has made effective field theory a very 
powerful tool for studying field theories at high temperatures. Perturbation theory breaks down for 
the resulting effective theory close to the phase transition and it is also severely infrared divergent 
in the symmetric phase. So one must use nonperturbative methods such as lattice simulations to 
investigate the phase transition. Dimensional reduction has been applied to a number of theories 



with spontaneously broken gauge theories: S77(2)+Higgs [27,29-30], the Standard Model |TL 



the MSSM [32-36], Sf/(5)+Higgs 0, the 2HDM @ and C/(1)+Higgs [27,38-39]. The three- 



dimensional effective theories have been studied numerically in e.g. Ref, [40-44] for S , f/(2)+Higgs, 
in Ref. |5| in the case of SU(2) x f/(l)+Higgs, in Refs. [39,46-53] for [7(1)+Higgs, and in Ref. §4 



for SU(3) x SU(2) with two scalar fields (The latter arises as an effective 3d theory for MSSM for 
some values of the parameters). 

In the present paper we reconsider the Two Higgs Doublet Model. The model is interesting 
in its own right, but the 2HDM (with the temporal component of the gauge field as an additional 



adjoint Higgs field) in 3d also arises as an effective theory for the MSSM [31]. In Ref. [|34]] di- 
mensional reduction for this model was carried out in the one-loop approximation. The second 
effective theory obtained by integrating out the timelike component of the gauge field was con- 
structed in Refs. |32|, |34l , also at the one-loop level. In order to obtain a satisfactory accuracy for 



the thermodynamic quantities describing the phase transition, the scalar mass parameters in the 



3d theory are needed to two-loop order [31|. This calculation is carried out in the present paper. 



The plan of the article is as follows. In section II we briefly discuss the Lagrangian for the 
Two Higgs Doublet Model. In section III we present the parameters of the first effective theory. 
Section IV is devoted to the scenario where one of the Higgs doublets is heavy, and is integrated 
out together with the timelike component of the gauge fields. In Section V we consider the case 
where both Higgs doublets are light and are retained in the final effective Lagrangian. Finally, 
in section VI we summarize. In Appendix A, the notation and conventions are given. We also 
list the necessary sum-integrals in the underlying theory as well as the three-dimensional integrals 
needed in the effective theory. In Appendix B, some details of the matching procedure are given 
by explicitly calculating a mass parameter in the first effective theory. 



3 



2 Two Higgs Doublet Model 



The Euclidean Lagrangian for the SU(2) gauge-invariant 2HDM without fermions reads 



+ ml($\$2 + $2*1) + V{$i, $2) • 



Here, $1 and $2 are the SU (2)-doublets 



!>2 + ir}2 



* 2 = 73 



!>3 + m 

i>4 + ^4 



and 



(d^-igT a A a j2) 



d,A a v - d v Al + 36*^ . 



(2) 
(3) 



Here, g is the gauge coupling, i = 1, 2 and r 1 , r 2 and r 3 are the three Pauli matrices. 
The potential V($i, $2) is 0, ^ Q 

y($!,$ 2 ) = A 1 ($ t 1 $ 1 ) 2 + A2($^2) 2 + A 3 ($ t 1 $l)(^$2) + A4($ t 1 $2)(^$l) 

+a 5 [($ t 1 $ 2 )($ t 1 $ 2 ) + + a 6 + ($ t 1 $ 1 )($ t 1 $ 2 ; 



+A 7 



(^$ 2 )($2$l) + ($^2)($ t l*2] 



(4) 



where the scalar self-couplings are denoted by Xi — X 7 . All calculations in the present paper are 
carried in the Landau gauge. This is merely a convenient choice, since many diagrams vanish 
in this gauge. We emphasize that the parameters of the effective Lagrangians are gauge fixing 
independent. 



3 Dimensional Reduction 

In this section we carry out the dimensional reduction step for the Two Higgs Doublet Model. 

The fields in the effective Lagrangian are identified (up to normalizations) with the zero- 
frequency modes of the fields in the full theory. If the fields in the effective theory are denoted by 
$' 2 , A\ and Aq, we can schematically write at leading order 

<^(x) « VT f dr$i(x, r) , A? (x) « VT f dr^ a (x, r) , <(x) ^ VT drA a (x, r) . (5) 
Jo Jo Jo 
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The effective Lagrangian consists of all terms which can built out of the fields $ 2 , A\ and Aq 
and which satisfy the symmetries present at high temperature. Examples of symmetries are three- 
dimensional gauge invariance and an 0(3) symmetry for the field Aq. The effective Lagrangian 
then reads 



cff 



\p'if*ij + (A^tA^) + (A* 2 ) f (A<f 2 ) + M?(aO*?*! + M 2 WM 

+ 1a £ (/,)(«) 2 + ^K/i)^^;^ + 4(m)$ 2 ^2« + < s - 



\af\2 



(6) 



Here, we have explicitly written the superrenormalizable part of the Lagrangian, while SC' eS rep- 
resents all higher order operators consistent with the symmetries of the Lagrangian. The gauge 
coupling is denoted by g|(/i), the quartic coupling constants in the potential V(<l >/ 1 ,$ 2 ) are de- 
noted by Aj(/i), and D^Aq = iPi + 9E^ abc A^)AQ. The parameters of C' eS encode the physics at the 
scale T and are called short- distance coefficients. The coefficients of the effective Lagrangian are 
determined by calculating static correlators in the full theory and calculating the corresponding 
correlators in the effective theory and require that they be equal at distances R^$> 1/T [27-28]. 

The matching procedure is complicated by the breakdown of the simple relations (|J). Beyond 
leading order we must allow for short- distance coefficients multiplying the fields in the effective 
theory. At the one- loop level (next-to-leading order), the short-distance coefficients are given by 
the momentum dependent part of the two-point functions, and are associated with field renor- 
malizations in the underlying theory. These parameters are called field normalization constants, 
and are denoted by SW'(O), ITgo and U^' for the Higgs fields, the timelike component of the 
gauge field and the spatial components of the gauge field, respectively. The relation between the 
fundamental scalar field $j and the scalar field $• in the effective theory can then schematically 
be written as 



£«'(()) 



1/2 



(7) 



and similarly for the other fields. The above remarks also apply when we consider the two effective 
three-dimensional field theories in the next section. 



Si(0) 



9g< 



64vr 2 



:L h 



nffi'(o) 



16tt 2 



4L h - 
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n«'(o) 



r 



16vr 2 



4A + 



The field normalization constants have been calculated and listed by Kajantie et al. in [31 
for the Standard Model with A" Higgs doublets. For N = 2, the results are 



(8) 



The coupling constants of the scalar fields have been calculated by Losada in |34| at the one-loop 
level. We list the results here for completeness. 

AiT - T [l2A? + XI + A 3 A 4 + \x\ + 2X\ + 6A* - \\ l9 * + ^] ^ + ^ 



Ai(/x) 



,(9) 
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A 2 (/i) 
As(/i) 

A 4 (/i) 

a 5 (aO 

A 6 (/i) 
A 7 (/x) 



A 2 T-T 
A 3 T-T 



12A 2 + A3 + A 3 A 4 + ^-A 2 + 2A 2 + 6A 2 - ^A 2 ^ 2 + ^ 
2 2 lb 



16vr 2 8 16tt 2 V J 



6A1A3 + 2AiA 4 + 6A 2 A 3 + 2A 2 A 4 + 2A 2 + A 2 + 4A 2 + 2A 2 + 8A 6 A 7 
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+ 2\ 2 7 --\ 3 g 2 + ~g 4 



L b + 3g 4 T 



A 4 T — T 

- 2^ 
A 5 T-T 



16tt 2 4 16vr 2 ' 
2AiA 4 + 2A 2 A 4 + 4A 3 A 4 + 2A 2 ; + 32A 2 , + 5\ 2 6 + 2A 6 A 7 + 5A 2 . 
U 

167T 2 ' 

4AiA 5 + 4A 2 A 5 + 8A3A5 + 12A 4 A 5 + 5A 2 ; + 2A 6 A 7 
U 



5A 2 7 



16tt 2 ' 

A 6 T - T [12AiA 6 + 3A 3 A 6 + 3A 3 A 7 + 4A 4 A 6 + 2A 4 A 7 + 10A 5 A 6 + 2A 5 A 7 
U 



9 x 2 



16tt 2 



A 7 T - T [12A 2 A 7 + 3A 3 A 6 + 3A 3 A 7 + 2A 4 A 6 + 4A 4 A 7 + 2A 5 A 6 + 10A 5 A 7 



9 x 2 



16tt 2 



(11) 
(12) 

(13) 
(14) 
(15) 



The coupling constants g%(/j), h 2 E (fi) and A^(/x) have been computed in e.g. Refs. |3T|, |34 

1 , 2 



h 2 E (ji) 
A S (/i) 



g 2 T 

-g 2 T 
4 y 

3g 4 T 
8tt 2 



1 + 



9 



16vr 2 



7U + - 



1 + 



9 



16tt 2 



7L h 



49 
~6~ 



3Ai A 3 A 4 



47T 2 47T 2 



(16) 
(17) 
(18) 



The coupling constants are all renormalization group invariant to this order, which can be verified 
by using the renormaliztion group equations. This property holds to all orders in perturbation 
theory if the effective Lagrangian is restricted to superrenormalizable terms. 

The scalar mass parameters have been computed in the one-loop approximation by Losada in 
Rcf. |M| . Here, we present results for the mass parameters to two-loop order: 



M 2 M 



9 



,2- 

rp2 



9 



6mf\i + Ira^X^ + m\\^ + Qm^\ — ^mlg 2 



16tt 2 



+ [6A 1 + 2As + A 4 +^]^ + 



T 2 r3, 2 1,2 1,2 45 4 
4 Xl9 + 4 W + 8 W+ 32* 



T 2 



167T 2 



12A 2 + 2A 2 + 2A 3 A 4 + 2A 2 + 12A 2 + 9A 2 + 3A 2 - 9A^ 2 - 3A 3 # 5 



6 



3 , 2 75 4 
2 Xi9 - 16 9 . 



, 3T 

m h c 



(19) 



M 2 (/i) 



m — 



9 



6/71^2 + 2m 2 A 3 + m 2 A4 + 6m 2 A7 — -g 2 ?™ 2 . 



16vr 2 



+ (6A 2 + 2A 3 + A4 + ^I)^+ 167r2 



T 2 T3 A 2 1 2 1 2 45 4 
^A 2 # 2 + -X 3 g 2 + -A 4 <? 2 + — £ 4 



T 2 



16tt 2 



12A 2 + 2A 2 + 2A 3 A 4 + 2A 2 + 12A 2 + 3A 2 + 9A 2 - 9A 2 # 2 - 3A 3 # 2 



-M 2 



75 
16*' 



, 3T 

m h c 

fi 



(20) 



M 2 (/i) 



m 3 - 



mgA 3 + 2mgA 4 + GmgAs + 3mf A 6 + Sm^Xj — ^g 2 m\ 



U 



16tt 2 



+ (A 6 + A 7 )- +16jr2 



rp2 

T 



T 2 rs 3 

J. O.o A 2 

o X a9 + „A 7 5f 



T 2 



167T 2 

9 



[6AiA 6 + 6A 2 A 7 + 3A 3 A 6 + 3A 3 A 7 + 3A 4 A 6 + 3A 4 A 7 + 6A 5 A 6 + 6A5A7 



9 



' 3T 
In h c 



Here, c is the constant 27 



9 C(2) 1 , 



-0.348725 . 



[22) 



Note that we have written our mass parameters in terms of the renormalization group invariant 
couplings of the 3d theory. The remaining dependence on fi reveals that scalar mass parame- 
ters depend explicitly on the scale fi. This dependence on the scale fi is canceled by the scale 
dependence arising from calculations in the effective theory. 



The Debye mass is normally needed in the one-loop approximation [31 

m|( M ) = g 2 T 2 . 
There is no dependence on fi at leading order in g 2 . 



(23) 



4 Integrating out Aq 

The next step is to integrate out the adjoint scalar triplet Aq. This is carried out by calculating 
correlators in the two theories at distances R ^> 1/gT and require that they be the same [27-28]. 
The parameters in the effective theory encode the physics on the scales T and gT and are called 
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middle-distance coefficients. Before doing this, however, we must determine the masses of the 
scalar doublets near the phase transition. This is done by constructing the scalar mass matrix 
and finding the temperatures at which it has zero eigenvalues. The higher of these temperatures 
is close to T c , where the phase transition takes place and the corresponding eigenvector (Higgs 
doublet) has a mass of order g 2 T. The mass of the second scalar multiplet (after diagonalization) 
is determined near T c and it is found that it is generally of order gT and it should be integrated 



out together with Aq [31]. Only with fine-tuning of the parameters in the 2HDM, is it possible to 



obtain a mass of order g 2 T [01] . In this case it must be kept in the second effective Lagrangian. 
Both cases are considered below. The diagonalization modifies the parameters of (|6|) and the 
relations between the old and new parameters can be found in [^, In the following it is the 
rotated parameters of (§) that appear in the formulas. 



4.1 One Heavy Higgs and one Light Higgs 



In this subsection we consider the case where one of the Higgs fields (denoted by $ 2 ) * s heavy and 
has a mass of order gT. Hence, we integrate out this field together with the adjoint Higgs field Aq. 
The second effective field theory is then SU (2) +one Higgs doublet with higher order operators 
satisfying the symmetries. The Lagrangian reads 



1 „~ 



4ff = ^G^. + M^)!)^ 



(24) 



The gauge coupling is denoted by g 2 M (n). 



The field renormalization constant for the scalar fields vanish in the one-loop approximation, 
since there are no momentum dependent one-loop diagrams with internal Aq"s. Thus 



(25) 



This is in contrast with the gauge fields, since there is a momentum dependent one-loop diagram 
with Aq on the internal lines and Af on the external lines. The result is [^ 



A? » A" 



1 



9% 



+ 



9% 



1/2 



2Anm E 48vrM 2 
The results for the coupling constants can be found in [32, 34 1: 

1 



(26) 



Ai(/i) = Ai(/i) 
9mM = £eO) 



3hl 



—- — [2kl + 2A3A4 + K + 4A^ + 24A^ - 24A 6 A 7 ] - 

lo7rM 2 QTrrriE 



1 - 



9 2 e 



9l 



24nm E 48ttM 2 



(27) 
(28) 



8 



The mass parameter has previously been computed by Losada [0 at one- loop. The result in the 
two-loop approximation is: 



M 2 M 



Ml(li)--r (2A 3 + A 4 ) M 2 + 3h 2 E m E 



47T 



+ 



9 



A 2 

1 



. u 1 
In — + - 

. Mo 2 



16tt 



-A 2 



In 



16n 2 

u 1 
1- - 

2. 



[2A 2 A 3 + A 2 A 4 + 6A 6 A 7 



3M 



16tt 2 



2A 2 + 2A 3 A 4 + 2A 2 + 12A 2 - 3A 3 g 2 E - ^A 4 <?|) In 



2M 2 



+ A 2 + A3A4 + At + 6A 2 - ^A 3 g 2 E - ^A 4 g 2 E 



16n 2 
1 

167T 2 



Qh% - \2h\g\ + -a% ) In + 3^ - 3^ 



,2 2 



6/l1 



M2 



2m£ 

,2 m E , 3. 2 m B 



+ 3A ^M7 + 2 A ^ 



Mo 



(29) 



4.2 Two Light Higgs Doublets 



In this subsection we consider the other scenario when both Higgs fields have masses of order g 2 T. 
The effective Lagrangian is now a three-dimensional 2HDM with additional higher order operators 
which satisfy the symmetries at high temperature: 



C cS = i<V4- + (A$i) f (A$i) + 

+ Ml(fi)(^ 2 + 14$!) + $ 2 ) + 5jC cS . 

The scalar couplings are denoted by A3 and the gauge coupling by g 2 M (n). 



(30) 



Again the scalar fields are not renormalized by integrating out the Aq fields and so fl25|) also 
holds in the present case. This is in contrast with the gauge fields, since there is a trilinear 
coupling between Aq and Af: 



A? « A" 



1 + 



9% 



-,1/2 



24irmE 

The scalar couplings Ai(/i)— A 3 (/i) get modified by the integrating out Aq: 



Ai(/x) = AM - 



A 2 (/i) = A 2 ( At ) - 



3^ 



AM = AM) - 



34 
47rmp 



(31) 



(32) 



The other coupling constants, A 4 (/i) — A 7 (/i), are not modified in this step. The gauge coupling 
reads 

9m{v) = 9e{v) 



9e 



247tm E 



(33) 
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The expression for the mass parameters at the two-loop level are 

3h E m E 1 



m 2 M 


= m 2 m 


M 2 2 (/x) 


= m 2 m 


mm 


= mm 



An 167T 2 
3h E m E 1 



47T 



16vr 2 



Qh E -12hW E + 3 -g E ) ln£- 



6/i 4 1 ^ 2 " 2 1 " 4 



+ 3/4 - 3/i|^| 



34 



,(34) 

,(35) 
(36) 



5 Summary 



In the present paper I have applied the effective field theory methods developed in Refs. [27-28] 
to the 2HDM. I have exploited the fact that there are three well separated momentum scales and 
constructed a sequence of two effective three-dimensional field theories. The parameters in the 
final effective Lagrangian have previosuly been calculated in the one-loop approximation The 
two-loop results presented here are new. 

The resulting field theory can be used for investigating several aspects of the phase transition 
in the Two Higgs Doublet Model. This includes in particular the strength of the phase transition, 
and also the sphaleron rate immediately after the completion of the phase transition. 



A Notation and Conventions 



Throughout the work we use the imaginary time formalism, where the four-momentum is P = 
(p , p) with P 2 = pi + p 2 . The Euclidean energy takes on discrete values, p = 2niiT for bosons. 
Dimensional regularization is used to regularize both infrared and ultraviolet divergences by work- 
ing in d = 4 — 2e dimensions, and we apply the MS renormalization scheme. We shall use the 
following notations for the sum-integrals that appear 



f p f{P) 



; 7 W f d 3 ~ 2e p 



47r / P0 =2*nT J ( 27r 



T E ll^JiP)- (a.i) 



The one-loop sum-integrals needed in this work have been calculated in e.g. Refs. [8,28,31]. 
We list them here for the convenience of the reader: 
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p (P 



2\2 



16tt 2 



- + L b + 0{e) 



p2 




P (P 2 ) S 



p2 



P (P 2 ) 5 



-^[l + e(/ -2)] 



+ L 6 + 2 + 0(e) 



P 1 



p (P 2 ) 4 



64tt 2 
1 

1287T 2 



Here 



2 In 



4vrT 



+ 2 7 i 



- + L h + - + 0(e) 



/ e = 21n^ + 2 7B -21n2 



'C(2) 



(A.3) 
(A.4) 
(A.5) 
(A.6) 

(A.7) 



Moreover, 7^ is the Euler-Mascharoni constant and £(x) is the Riemann Zeta function. 

In the calculation of the mass parameters, we also need the fact that two-loop setting sun 
diagram is zero |2"8fl: 

' . (A.8) 



fpQ P 2 Q 2 (P + Q) 2 

In the effective three-dimensional theory we also use dimensional regularization in 3— 2e dimensions 
to regularize infrared and ultraviolet divergences. In analogy with ( |A.1| ), we define 



4tt J 



d 3-2e 



P 



(2< 



3-2e 



HP) 



(A.9) 



Again \i coincides with the renormalization scale in the modified minimal subtraction renormal- 
ization scheme. 



The one-loop and two-loop integrals needed are 

1 



p p z + m z 



m 

47T 



[1 + 0(e)] , 



1 



P (p 2 + m 2 )(p 2 + m 2 ) 
1 



An {mi + m 2) 



[1 + 0(e)] 



/ 

J P q (p 2 + m 2 )(q 2 + m|)[(p — q) 2 + rn 2 



1 



167T 2 

+ ln 



1 1 



mi + m,2 + ms 



+ 0(e) 



pg (p 2 + m 2 )(q 2 + m 2 ) 2 (p — q) 2 



l&Ti 2 m 2 
1 



pq (p 2 + m 2 )(q 2 + m 2 )(p — q) 



167r 2 m 2 



4 + ^ 
4 + 



(A.10) 
(A.ll) 

(A.12) 
(A.13) 
(A. 14) 



These integrals have been computed by several authors, e.g. in Refs. [pL |28|, ^1 
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B Matching Example 



In this appendix we explicitly show how the matching procedure is carried out by determining the 
mass parameter M 2 (A) to two- loop order. 

We denote the static two-point function of the Higgs field in the full theory by r^^k), and 

(2) 

the static two-point function in the effective theory by 1^, ^ (k). The corresponding self-energies 

are denoted by E(k) and S(k). Finally, the n'th order contribution to the self-energies in the loop 
expansion are denoted by £( n )(k) and E( n >(k). 

The self-energies can be expanded in powers of the external momentum k and so we can write 
the two-point functions as 



r^(k) = fc 2 + m2 + S«(0) + A; 2 S( 1 >(0) + S( 2 )(0), 



rg^(k) = fc 2 +M 1 2 (^)+E( 1 )(0)+A ; 2 E( 1 ) / (0)+E( 2 )(0)+5M- 



(B.l) 

:b.2) 



Here, we have added a mass counterterm, which is associated with mass renormalization. The 
mass parameter is then determined by matching these two-point functions, and by taking the field 
normalization constant into account, we can write the matching equation as 



r£ ) ^(k) = [i + E«'(o)] rg^(k). 



(2) 



(B-3) 



Since the external momentum k provides the only mass scale in the loop integrals contributing to 
the self-energy of the effective theory, they all vanish in dimensional regularization. The matching 
equation ( |B.3| ) can then be rewritten as 

M 2 Gu) = m\[l- E (1) '(0)] + E (1) (0) [l - E (1) '(0)] + E (2) (0) - 5M 2 . (B.4) 
The self-energy at one-loop order in the full theory reads: 



E«(k) 



1 

p P 1 



6m 2 Ai + 2m 2 A 3 + m 2 A 4 + 6m 2 A 
+ [6Ax + 2A 3 + A 4 + - A {d - 1)/] ^ ^ 

_ 3o ^ k l + 3a *4 (P k ) 2 

y Tp P 2 (P + Kf y Tp P\P + Kf 



(B.5) 



Expanding in powers of the external momentum k gives 



E«(k) 



+ 



6771^1 + 2m 2 A 3 + m 2 A 4 + 6m 3 A 6 | ^ — 



6A X + 2A 3 + A 4 + \{d - 1)/] f p ±- ~/k 2 f p ± + 0(k'/T 2 ) . (B.6) 
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This implies 



S (1) (0) 



6m 2 Ai + 2m 2 ,X 3 + m^X^ + 6m 2 A6 



1 

p P 1 



[6A 1 + 2A 3 + A 4 + ^(rf-% 2 ]^-^ 
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73 d 



1 



4" fpP 4 ' 

The two-loop contribution to the self-energy at zero external momentum is 



(B.7) 
(B.8) 



£( 2 )(0) 



36A 2 + 12AiA 3 + 6A1A4 + 12A 2 A 3 + 6A 2 A 4 + 4:\j + 4A 3 A 4 + X\ 
+ 18A 2 + 18A 6 A 7 + ®-{d - l)X l9 2 + ^(d- l)X 3 g 2 + ^(d- l)X 4 g 2 



+ l(^-9d + 6) 9 ^ pQ ^ QA . 



(B.9) 



After renormalization of the mass parameter m\ as well as the coupling constants, we are left with 
a pole in e. This pole is canceled by the mass renormalization counterterm, which is 



5Ml 



12X\ + 2A 2 , + 2A3A4 + 2A^ + 12A 2 , + 9A 2 , + 3A 2 . - 9X l9 2 - 3A 3 ^ 2 
_ 3 2 _ 75 4 1 1 
2 49 16 9 I 647r 2 e ' 



(B.10) 



This is the result for the mass counterterm of the three-dimensional 2HDM at next-to-leading 
order in the coupling constants. The mass parameter Mf(fi) is then given by (|19|). 

This work was supported in part by a Faculty Development Grant from the Physics Depart- 
ment of the Ohio State University and by the Norwegian Research Council (project 124282/410). 
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